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Abstract: - Wavelet collocation method for numerical solution nth order Volterra integro diferential equations 
(VIDE) by expanding the unknown functions, as series in terms of chebyshev wavelets second kind with 
unknown coefficients. The aim of this paper is to state and prove the uniform convergence theorem and 
accuracy estimation for series above. Finally, some illustrative examples are given to demonstrate the validity 
and applicability of the proposed method. 
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I. INTRODUCTION 

Basic wavelet theory is a natural topic. By name, wavelets date back only to the 1980s, on the 
boundary between mathematics and engineering, wavelet theory shows students that mathematics research is 
still thriving, with important appUcations in areas such as image compression and the numerical solution of 
differential equations [1], integral equation [2], and integro differential equations [3,8]. The author believes that 
the essentials of wavelet theory are sufficiently elementary to be taught successfully to advanced 
undergraduates [4]. 

Orthgonal functions and polynomials series have received considerable attention in dealing with 
various problems wavelets permit the accurate representation of a functions and operators. Special attention has 
been given to apphcation of the Legendre wavelets [5], Harr wavelets [6] and Sine-Cosine wavelets [7]. 

The solution of integro-differential equations have a major role in the fields of science and ingineering 
when aphysical system is modeled under the differential sense, it finally gives a differential equation, an 
integral equation or an integro-differential equations mostly appear in the last equation [8]. 

In this paper the operational matrix of integration for Hermite wavelets is derived and used it for 
obtaining approximate solution of the following nth order VIDE. 

u<"\x)=g(x)+/; k{x, t)u^'\t)dt (1) 
where n > s k(x,t) and g(x) are known functions, and u(x) is an unknown function. 

II.SOME PROPERTIES OF SECOND CHEBYSHEV WAVELETS 

Wavelets constitute a family of functions consti'ucted from dilation and ti-aslation of a single function 
(»(t) called the mother wavelet. When the dilation parameter a and the translation parameter b vary continuously 
we have the following family of continuous wavelets as [9]. 

<»a,i(t) = |a|l"o(^),a,b £ R, a ?i 0. 
The second chebyshev wavelets ffl^„,m(t) = (a(k,n,m,t) involve four arguments, 
n = 1, ... , 2*^"^ , k is assumed any positive integer, m is the degree of the second chebyshev polynomials and t 
is the normahzed time. They are defined on the interval [0,l)as 

<.(t) = f2^ f^™(2'^t-2n + l),^<t<^ (2) 

^ 0 otherwise 
where U^{t) = Jl^mCO m = 0,1, ... , M - 1 (3) 

here U^(t) are the second chebyshev polynomialsof degree m with respect to the weight function w(t) = 
Vl — on the interval [-1,1] and satisfy the following recursive formula 
f/o(t) = l,t/i(t) = 2t, 
U^+iit) = ZtU^it) - U^_iitl m = 1,2,... 

The set of chebyshev wavelets are an orthonirmal set with respect to the weight function 
M/„(t) = w„(2''t-2n-l-l) . 
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III.FUNCTION APPROXIMAT 

A function /(t) defined over [0, 1) may be expanded as 

/(t) = K=lI:S=oC„n^«^n^(t) (4) 

Where = (/(t), oi^ it)) 

In which (. ,.) denoted the inner product in L^^[0,1). If the infinite series equation (2.34) is truncated, 
then it can be written 

fit) = iti l^-=l C^m dm(t) = C^a/ (t) ^ (5) 

C = [Qo'Qi' —>^1{M-1)> ^20> — >^2iM-l)> ■■■,^2^-1, ... ,C2k-l,^_i] fi^ (t) = 

[4o (t), fi^i (t) (t), (4o (t) (t). - fi|.-lo (t). - ft^z^-lM-l )]^ 

Convergence Analysis 

for Chebyshev wavelets of Second kind. 

A function f(t) defined over [0,1] may be expanded as: 

m = Y.n=li:^=ofnmC^nmit) (6) 

where 

fnm = im,aimit)) (7) 

In eq.(5).(. , . ) denotes the inner product with weight function m/„ (t). 
If the infinite series in eq.(4) is trancated then eq.(4) can be written as: 

/a)=/2^-lM-l W = n = l YIt=lfnm 

where F and ai-e 2^M Xl matrices given by 

P = [/lO-Al fmJlO /2M-1 /z^o fi^-^M-iY (8) 

o/(t)=[4o (t), ft^i (t) 4m-i (t), c4o (t) c4m-i (t), W '"a^-iM-i 

Theorem(l): (Convergence Analysis theorm) 

Assume that a function f(t) £ L^. [0,1], w* = Vl — with |/"(t)| < L, can be expanded as infinite series 
of second kind chebyshev wavelets, then the series converges uniformly to f(t). 
Proof: since = (/(t), (0) 
then 



Jo 



fnm = nt)oi^{t)w^mt 



=L-i, ' ^/(t)f^m(2''t-2n + l)w(2''t-2n+l)dt 

/2k-l 

If we make use of the substitution Z'' t — 2n + 1 = cos 6 in (5), yields 
V2 



V2 r'' /cos0 + 2n-l\ 

/nm=n^ / sin0sin(m+l)0d0 

22V^^o ^ ^ ^ 



By using the integration by pai'ts. 
Then eq.(lO) becomes 

(COS e + 2n-l\ /sin mO sin(m + 2)0 V 



1 [ /COS0 + 2n - 1\ /si 



1 C /cos0 + 2n-l\ 

=--^ /' sin.sinmed. 

mv2 2 2 Vtt 0 



1 /cos0 + 2n-l\ /; 



+ 2 

sinmfl sin(7n + 2)0\ 



1 /cos0 + 2n-l\ 

— /' 1 sin0sin(7n + 2)0d0 

(m + 2)V2 2XV^^o V 2 / 



After performing integration by parts again yields. 
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^ /COS0 + 2n — 1\ TcosCm — 2)0 — cosmfl cosmS — cos(m + 2)0l 



m — 1 m + 1 



1 /COS0 + 2n - 1\ rcos(7n)0 - cos(m + 2)0 

3 5fc I / ( Tfc ) ^ . 1 

(m + 2)22 2XV^ ^0 ^ ^ n "^+^ 



(m + 2)22 2XV^ ■'o 
cos(7n + 2)9 - cosim + 4)0] 



If" /COS0 + 2n - 1\ [( 



(12) 

cos(7n - 2)0 - cosm0 cosmS - cos(m + 2)0l 1^ 
^1 m + 1 J I 

cos(m - 2)0 - cos m9 cos m0 - cos(7n + 2)0] I 



If'' /COS0 + 2n - 1\ [cc 
" f k" C^°^^ ~ ''"^|^d0 X J" 1*''^ ~ '^"^'■"^ - 2mcosm0 + (m + 1) cos(m - 2)0 1^ 



Jo 



(m- l)(7n + 1) 

"(m- ly cos^Cm + 2)0 + 4m'^ cos^ mO + {m+ lY cos^{m - 1) 
(m-l)2(m + l)2 

im-mm + iy k - 1) + + 2 + 1) 1 

[3m2 + 1] 



(m-iyim + iy 

Thus, we get 

/cos 0 + 2n - 1\ [cosCm - 2)0 - cos ?n0 cos m0 - cos(7n + 2)0 



If'' /COS0 + 2n- 1\ [( 



m-1 m+1 
nL(_3m^ + lY^^ 
(m- l)(m+ 1) 

- 1\ rcos(7n)0 - cos(7n + 2)0 cos(7n + 2)0 - cos(7n + 4)0] 



(13) 



+ 3 

'COS0 + 2n - 1\ 



2" J 

(m + 1) cos(m + 4)0 - (m + 1) cos(7n + 2) 0 + (m + 3) cos(m + 2)0 + (m + 3) cos(m)0|^ 
^ (m + l)(m + 3) I 

fl /cos0 + 2n-lM^ 

4. r( — 3 — 

■ |7ncos(7n + 3)0 - i2m - 2) cos(7n + 1)0 + (m + 20)cos(7n - 1)0 1^ 



m{m + 2) 

f " (m + 1)^ cos2(m + 4)0 + (2m + 4)^ cos2(7n + 2)0 + (m + 3)^ cos^ (771)0 
(m+ l)2(m + 3)2 

[6m^ + 24m + 26] 



2(m + l)2(m + 3): 

= 7 TT^JT ^ [Sjti + 12m + 13] 

(m+ l)2(m + 3)2 

Thus we get 

■ f " „ /cos 0 + 2n - 1\ rcos(m)0 - cos(m + 2)0 cos(m + 2)0 - cos(m + 4)0] | 



If" /COS0 + 2n- 1\ [cc 



m + 1 m + 3 

-(3m2 + 12m + 13)^/2 



(m + l)(m + 3) ^' 

(14) 

Using eqs.(13)and(14), one can get 
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5k-3 ( nLi3m^ + 1) ^2 nL{3m^ + 12m + 13) 
\fn^\<2 2 n 2 \^ _ ^^^^ ^ ^) - ^ 2)(m + l)(m + 3) j 



22(m+l) ^^"""^^ ^ 22(m-l)2 

Finally since n < Z'' — l,then 



Accuracy Estimation of fi^^ W= 

If the ftinction f(x) is expanded interms of fourth kind chebyshev wavelets, 

/(X)= K=li:^=0 /nmfi^^W (15) 

It is not possible to perform computation an infinite number of terms, therfore we must truncate the series 
in (15). In place of (5), we take 

fm W = Y.t=l YJt=l fnm dm W (16) 

SO that 

fix) = (X) + I.m=M fnm dm C^) 

or /(x)-/M(x) = r(x) (17) 

where r(x) is the residual function 

r(x) = i::=2^-l+l ^m=M fnm dm (^) (18) 

we must select coefficients in eqs.(17) and (18) such that the norm of the residual function ||r(x)|| is less than 
some convergence criterion £,that is 

for all M greater than some value Mg. 
Theorem (2) 

Let f(x) be a continuaus function defined on [0,1), and |/ (x)| < L, then we have the following accuracy 
estimation 

c.,.<^i::_,.-.,,i:s..^-^ (19) 



where 



Proof:- 



c,,M = {f,\rix-)fwAx-)dxy 



Since Cfc^ = [f^{rix)fw^ix) dxf 

CL=f,\rix)fwAx-)dx 

=fo I]^=2''-l+l ^m=M fnm (o/w^ (x))dx =I]^=2''-1+1 ^m=M fnm Jq (dm (^))^ W„ (x) dx 

or 

CL=I.:=2^-l+lI.m=M JoV4(fi^)'w„(x)dx 

CL = I.n=2k-i+iI.m=M f^ (z^) 0^(2" X - 2n + l)Vl- (2'^x-2n+ l)2dx 

Lett = 2''x-2n + l, 

we have. 



(t)Vl-t^dx = 



Cfcm — I]„=2''-l+l ^m=M f^ 2 
Clem < I]^=2''-l+l^m=M 8(m-l)4 
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wavelet collocation method for VIDE with mth order: 

In this section the introduced wavelets collocation will be applied to solve VIDE with mth order, 
uf^ W = 91 W + Jo" K^j (x, t)uf it)dt ,n>s (20) 

With the following conditions uf (0) = a^^ i = 1,2, ... ,1 s = 0,1,2, ... ,n — 1 

AfUnction uf (x) which is defined on the interval x £ [0,1] can be expanded into the second chebyshev 

wavelet series 

<(x) = i:f=iqfi;i(t) (21) 
Where C; are the wavelet coefficients. 
Integrate eq.(21) m times,yields 

u(x) = Sflo ct Jo' ... Jo' CO, (t)dt + YT=~o' j; dm-i (22) 
Using the following formula 

f ... f (o,{t)dt = f (x - t)"-'®i(t)dt 

therefore eq.(22) becomes 

u(x) = i:f=o Ci ^ Jo'(x - t)"-^ »i {t)dt + 1:7=0' 7^ (23) 
Let (x, t) = ^^^rj^ and = Jj" /C„ (x, 0^; (t)dt i=0, 1, . . .,M 
This leads to 

In similar way, we can get 

u(^) (x) = Xflo c^ir' + 1:7=0'-' (24) 
Substituting eqs (22) and (24) in (20), yield 

Sf=i Q (t) = 5; W + /; /f;,- (x, t) [Sflo qZ-r^ + i:7=o'-i a„_,_^. ] dt (25) 
or Ylt=i Ci (t) - (^) = 5; (^) + T.'IZt^ ^ (^) (26) 

where A^{x) = Kn{x,t)Ul-' {t)dt i=0,l,2,...,M 

S. (x) = /J" /C„ (x, t) t-* dt j=0, 1 ,2, . . . ,n-s- 1 (27) 

Next the interval x £ [0,1] is devided in to I Ax = y and introduce the collocation points 
Xfc = '^-j- , k=l,2,...,l eq(21) is satisfied only at the collocation points we get asystem of linear equations 
i:r=i Ci [ffiS (X) - A; (x)] = gi (x) + Sp-' ^ Bj (X) (28) 

The matrix form of this system 
is C F=G+ I]f=o ^^y^ Bj (x) where F=o (x), G=g(x) 

1. Design of the matrix A:- 

When chebyshev wavelets second kind are integrated m times, the following integral must be 

evaluated. 

VI = (x, t)®; it)dt , i=0, 1, 2, . . ., M 



" 2*(n-l)! 



0 0 0 0 
0 0 0 0 



Therefore the matrix Aj (x) can be constructed as follows 
Since (x) = Jj" /C„ (x, t)^-' {t)dt i=0, 1,2,. . .,M 
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I K^{xo,t)Ul-'{t)dt 1 = 0 
Jo 

I K^{xi,t)Vl-' {t)dt i>0 

--'0 



IV. WAVELETS METHOD FOR VIDE WITH NTH ORDER 

For solving VIDE with mth order the matrix L" (x) in section above will be followed to get 
SHi Ci [o\ (xl - AJ] = gCxJ + Sp-i ^ Bj (xj I 6 [a, b] 
But^i(Xi) = JJ''-/C„(xL,t)i'?"^(t)dt where i=0,...,M 
Bj(xL)=Ci^„(xL,t)t"-^dt 
where Z,""^(t) as in eq(17),(18) 
that is Ai (Xi) = ^z, , F; (x^) = (x^) - (x^) = 
Numerical Results: 

In this section VIDE is considered and solved by the introduced method, parameters k and M are 
considered to be 1 and 3 respectively. 

Example 1: Consider the following VIDE: 

U"(x) = e2--j;e2(-t)t/'(t)dt 
Initial conditions U(0) = 0,U '(0) = 0. 

The exact solution f/(x) = xe^ — + 1. Table 1 shows the numerical results for this example with k=2, 
M=3 with error =10"^ and k=2, M=4, with error =10"'' are compared with exact solution graphically in fig. 



Table 1 : some numerical results for example 1 



X 


Exact solution 


Approximat solution 


Approximat solution 






k=2,M=3 


k=2,M=4 


0 


0.00000000 


0.00000001 


0.00000001 


0.2 


0.02287779 


0.02280000 


0.02287000 


0.4 


0.10940518 


0.10945544 


0.10940544 


0.6 


0.27115248 


0.25826756 


0.27826756 


0.8 


0.55489181 


0.54330957 


0.55330957 


1 


1.00000000 


0.99999995 


0.99999998 




Fig 1: Approximate solution for example 1 

Example 2: Consider the following VIDE : 

U(5)(x) =-2sinx + 2cosx-x + j^ix - f)U^^\t)dt 
Initial conditions t/(0) = 1, f/'(0) = 0, f/"(0) = -1, U^iO) = 0. 

The exact solution U (x) = cos x. Table 2 shows the numerical results for this example with k=2, M=3 with 
error=10"'' and k=2, M=4, with error =10"'* are compared with exact solution graphically in fig, 2. 
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Table 2: some numerical results for example 2 



X 


Exact solution 


Approximat solution 
k=2,M=3 


Approximat solution 
k=2,M=4 


0 


1.00000000 


0.99812235 


0.99999875 


0.2 


0.98006658 


0.98024711 


0.98005541 


0.4 


0.92106099 


0.92158990 


0.92104326 


0.6 


0.82533561 


0.82479820 


0.82535367 


0.8 


0.69670671 


0.69689632 


0.69678976 


1 


0.54030231 


0.54032879 


0.54035879 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Fig 2: Approximate solution for example 

V. CONCLUSION 

This work proposes a powerful technique for solving VIDE second kind using wavelet in collocation 
method comparison of the approximate solutions and the exact solutions shows that the proposed method is 
more faster algorithms than ordinary ones. The convergence and accuracy estimation of this method was 
examined for several numerical examples. 
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